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We have used first-principles density-functional-theory methods with a random-structure- 
searching technique to determine the structure of the previously unidentified post-cotunnite phase 
of Te02. Our calculations indicate a transition from the cotunnite to post-cotunnite phase at 130 
GPa. The predicted post-cotunnite structure has P2\/m space group symmetry and its calculated 
x-ray difi^raction pattern is in good agreement with the available experimental data. We find that 
the cotunnite phase re-enters at about 260 GPa. 

PACS numbers: 62.50.-p,71.15.Nc,61.50.-f,91.60.Hg 
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INTRODUCTION 

The majority of matter in the solar system is subject 
to pressures above 10 GPa,fli which motivates studies of 
materials such as oxides at high pressures. In the case of 
AX2 compounds, where A and X are, respectively, a di- 
valent cation and halogen atom or a tetravalent cation 
and oxygen atom, the general effect of increasing the 
pressure is to distort the anion polyhedra and eventually 
to increase the coordination number (CN) . [2] The highest 
CN observed in metal dioxides is in the PbCl2-type co- 
tunnite structure with CN=9. Metal dioxides with large 
cation radii often form cotunnite phases at high pres- 
sures, such as TiOi, Zr02, Hf02, Ce02, Pb02, PUO2, 
UO2, Tb02, Te02,j3i and Th02.[i The very important 
oxide Si02 is predicted to adopt the cotunnite structure 
above 690 GPa, which may be relevant to the study of 
extrasolar planets. fB^I The hardest known oxide is the co- 
tunnite structure of Ti02, which has been synthesised at 
high pressures and recovered to ambient conditions. @ 

Materials that adopt the cotunnite structure are ex- 
pected to transform under additional applied pressure 
into post-cotunnite structures with an accompanying 
increase in CN to 10 or more, as reported for some 
dihalides.Q In reviewing the high pressure phases of 
dioxides our attention was drawn to Te02 which, to the 
best of our knowledge, is the only dioxide for which a 
transition to a post-cotunnite phase has been observed, [sl 
Sato et a/. [8] studied Te02 up to pressures of 150 GPa in a 
diamond anvil cell. X-ray diffraction data showed strong 
evidence for a structural phase transition around 80-100 
GPa, but the quality of the data was insufficient to allow 
a determination of the structure of the new phase, al- 
though the known post-cotunnite structures of dihalides 
were eliminated. [8] Identifying the post-cotunnite struc- 
ture of Te02 would further our understanding of dioxides 
at high pressures. 



Te02 is also an interesting material from the point 
of view of fundamental science and technology. ,9] It has 
shown promise as a material for nonlinear optical devices, 
usually in a glassy form but potentially from nanosize 
crystals. [13, [nl 



RANDOM STRUCTURE SEARCHING 

First-principles ov Ah Initio (AI) Density- Functional- 
Theory (DFT) methods have been widely applied to ma- 
terials at high pressures, and have provided both confir- 
mation of experimental results and predictions of new 
phases and their properties. DFT calculations have 
given very accurate descriptions of the high-pressure 
phases of sp bonded materials. [l3| We have studied high- 
pressure phases of Te02 using AI DFT methods com- 
bined with "Random Structure Searching" (the AIRSS 
approach) . [isj This approach has been used to pre- 
dict high pressure phases which have subsequently been 
found experimentally, [13, 14 1 and to predict new high- 
pressure phases of materials such as hydrogen [l3 and 
ammonia. 



The recipe for a random search commences by generat- 
ing a set of initial structures, for each of which a random 
unit cell is created and renormalised to a reasonable vol- 
ume, and the desired number of each atomic species is 
randomly distributed throughout. Each of these initial 
configurations is relaxed to a minimum in the enthalpy at 
a predefined pressure, and the procedure repeated until 
the lowest enthalpy structures have been found several 
times. Such random searching is largely unbiased, but it 
can often be made much more efficient by applying con- 
straints. Any reasonable structure of Te02 will contain 
only Te-0 bonds, and therefore we have performed most 
of the searches by placing 0-Te-O molecules within the 
cells, rather than separate atoms. Another constraint we 



2 



1^] 


U. 1 o 






o 


I' ' ' 




0) 
H 






P2j /m 


per 




\\ 
\* 


pj 


U. 1 U 


Lmcm 


> 








03, 




w 




rence 




»\ 
\\ 


' • • , 


CD 




\\ 






0.00 


\\ 




alpy d 


\\ 






-0.05 






LU 









100 150 200 250 300 350 400 
Pressure (GPa) 

FIG. 1: (Color online) Enthalpy per Te02 unit relative to 
that of the cotunnite structure, as a function of pressure. 



have employed is to reject initial configurations in which 
atoms are closer than a defined minimum separation. We 
have also generated initial configurations with the space 
groups which contain a specified number of symmetry 
operations, and then relaxed the structures while main- 
taining the symmetry. 

Our DFT calculations were performed using the 
CASTEP 17] plane wave code with the Perdew-Burke- 
Ernzerhof (PBE) generalised gradient approximation 
(GGA) exchange-correlation functional and ultrasoft 
pseudopotentials.[3] For the searches we used a plane 
wave cut off energy of 490 eV and a Monkhorst-Packfioj 
Brillouin Zone sampling grid spacing of 27r x 0.07 A^^. 
All of the results reported in this paper were obtained by 
refining the structures obtained in the searches and cal- 
culating their properties using a higher level of accuracy 
consisting of a plane wave cut off energy of 800 eV and 
a grid spacing of 27r x 0.03 A~^. 

We first performed searches at 150 GPa. Uncon- 
strained searches were performed using 2 and 4 formula 
units of Te02. Another set of searches was performed us- 
ing initial configurations built by applying the symmetry 
operations of space groups chosen randomly from those 
with n operations to the randomly chosen positions of 
a Te atom and two O atoms, with n = 3, 4, 6 and 8, 
all subject to a minimum separation of rmin = 1.3 A. 
Searches were then performed using 0-Te-O molecules 
with initial bond angles of 120°, with 1, 2, and 3 molec- 
ular units and space groups with n = 4 operations, again 
with Tniin — 1.3 A. Additional searches were performed 
at 280 GPa using molecules with initial bond angles of 
120°. We used 1 molecular unit with n = 4 symmetry 
operations, 3 molecular units and n = 2 symmetry opera- 
tions and Tmin — 1-3 A, and 4 molecular units with n = 1 
symmetry operations and Vyam = 1-2 A. The searches 
produced a total of about 1800 relaxed structures. 
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FIG. 2: (Color online) Comparison of observed (red solid 
lines) and calculated (green dashed lines) x-ray diffraction 
data for the cotunnite structure (upper box) at 53 GPa, and 
the P2i/m structure (lower box) at 130 GPa. The exper- 
imental diffraction data is from Ref. [^. The experimental 
and calculated data were obtained with an x-ray wavelength of 
A = 0.4254 A. The black arrows mark known impurity lines in 
the experimental data,g] although other impurity lines may 
also be present. 



RESULTS FROM STRUCTURE SEARCHING 

Enthalpy-pressure curves for the more stable phases 
are shown in Fig. [TJ A structure of P2i/m symmetry 
was consistently the lowest enthalpy phase found at 150 
GPa in all searches with a total of 4 formula units, and 
also in the 8 unit search with 2 symmetry operations ap- 
plied to 4 molecular units. The structure with the second 
lowest enthalpy in these searches was always found to be 
the Pnma cotunnite structure. Searches with 2 and 3 
formula units did not yield structures with enthalpies as 
low as the searches with 4 or more, but a search with 
6 formula units produced a low-symmetry PI structure 
which has an enthalpy between that of cotunnite and 
P2i/m at 150 GPa. The cotunnite, F2i/m, PI and 
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FIG. 3: (Color online) a) Cotunnite, b) P2i/m, c) Cmcm and 
d) PI structures. The yellow (light) spheres are Te atoms and 
the red (dark) spheres are O atoms. 

Cmcm structures are shown in Fig. [3] and their struc- 
tural parameters are reported in Table HI The 12 unit 
search did not reveal any new structures that were lower 
in enthalpy than those already mentioned, although a 
previously unseen and fairly-low-enthalpy structure with 
space group P2i/c was found. The searches at 280 GPa 
did not yield any further low-enthalpy structures. 

In our experience, the appearance of a very low sym- 



metry structure, such as PI, as a low enthalpy phase 
suggests that another structure of even lower enthalpy 
might exist. We therefore performed an additional type 
of search using the cell obtained by doubling that of PI in 
each direction, giving a cell containing 48 formula units. 
We then performed "shakes" of the larger structure in 
which all atoms were displaced in random directions by 
a distance chosen randomly between and 0.25 A, and 
then relaxed, but in each case the original PI structure 
was recovered. 

The transition from cotunnite to P2i/m occurs at 130 
GPa in our calculations. Sato et al. observed a phase 
transition at 80 GPa after heating the sample to 1000 
K, and at 100 GPa at room temperature. Heating helps 
in overcoming kinetic barriers which are expected to be 
large in oxides and can also help to reduce anisotropic 
stresses. A temperature of 1000 K could, however, af- 
fect the coexistence pressure. The agreement between 
the measured transition pressure and the theoretical co- 
existence pressure is satisfactory, given the uncertainty 
in the experimental transition pressure and the fact that 
our calculations are at zero temperature. The maximum 
stabilization of the Plxjm phase over cotunnite is about 
0.031 eV per Te02 unit at 175 GPa, which is quite small 
but easily resolved in our calculations. Such small en- 
thalpy differences are often given quite accurately in DFT 
calculations for sp bonded materials where the volumes 
and the nature of the inter-atomic bonding in the two 
phases are very similar, as is the case here. 

Fig. [1] shows that the PI structure is marginally the 
most stable in the pressure range 248-269 GPa, although 
the enthalpies of the PI, P2i/m, and cotunnite phases 
differ by less than 0.0023 eV per formula unit in this 
range. The cotunnite structure becomes more stable 
than the PI and Plxjm structures again at around 260 
GPa. This re-entrant behavior of the cotunnite struc- 
ture is quite unexpected. The origin of the apparent 
'kink' in the enthalpies of the other structures relative 
to cotunnite in Fig. [T] actually lies with the nature of the 
cotunnite structure itself, as highlighted in Fig. [D At 
pressures up to about 160 GPa, the compressibihty of 
the cotunnite phase is nearly constant and is larger than 
that of P2i/m, but at higher pressures the compressibil- 
ities are similar. The region of high compressibility of 
the cotunnite structure is predominantly associated with 
compression along the a axis, whilst the c axis actually 
increases in length from 125 to 150 GPa, before contin- 
uing to decrease steadily with increasing pressure. The 
cotunnite phase has a larger volume at pressures below 
about 160 GPa, but it is slightly smaller at higher pres- 
sures, which tends to favor cotunnite over P2i/m. 

The theoretical and experimental diffraction data for 
the cotunnite structure shown in Fig. [2] are in very good 
agreement. The discrepancies in relative peak heights 
might arise from the form factors used to generate the 
theoretical data, from the differences in structures due 
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Space group 




Lattice parameters 

(A,°) 






Atomic coordinates 
(fractional) 




Pnma 


a=4.927 


6=3.223 


c=6.389 


Tel 


2398 


9500 


0.6104 




a=90.00 


/3=90.00 


7=90.00 


01 






0.9378 










02 


0.0482 




0.3035 


P2i/m 


a=6.287 


&=3.577 


c=4.475 


Tel 


0.1131 


0.7500 


0.2227 




a=90.00 


/3=97.15 


7=90.00 


Te2 


U .OU\JO 




0.7361 










Ol 


0.0711 


0.7500 


0.6543 










02 


2fi38 


9500 


0.2032 










03 


3930 


0.7500 


0.4901 










04 


0.3821 


0.7500 


0.9774 


Cmcm 


0=3.014 


6=10.144 


c=3.232 


Tel 


1.0000 


SSOfi 


0.7500 




a=90.00 


/3=90.00 


7=90.00 


01 


1.0000 


0.7527 


0.2500 










02 


1.0000 


0.5783 


0.7500 


PI 


a=4.473 


6=5.963 


c=6.280 


Tel 


0.1777 


8399 


0.1146 




a=99.95 


/3=97.56 


7=111.84 


Te2 


0.5173 


0.4923 


0.2671 










Te3 


0.8648 


0.1658 


0.4231 










Ol 


0.1107 


0.1638 


0.1524 










02 


0.0904 


0.5140 


0.3393 










03 


0.4251 


0.1584 


0.4611 










04 


0.6011 


0.1671 


0.1441 










05 


0.7083 


0.8322 


0.1993 










06 


0.7564 


0.5003 


0.0060 



TABLE L Structures of the cotunnite {Pnma, Z = 4 formula units per primitive cell), P2\/m (Z = 4 formula units per 
primitive cell), Pi {Z — 6 formula units per primitive cell), and Cmcm (Z = 6) phases of Te02 at 130 GPa. 
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FIG. 4: (Color online) Volume per Te02 formula unit of the 
cotunnite and P2i/m structures. 



to the approximate DFT and from the lack of temper- 
ature effects in the theoretical structure. The black ar- 
row indicates an impurity line in the experimental data 
identified by Sato et a/.,[8| which we see is absent in the 
theoretical data. The level of agreement between the 
theoretical and experimental x-ray data for the cotun- 



nite phase gives us a benchmark for making a similar 
comparison for the post-cotunnite phase. The good level 
of agreement between the theoretical and experimental 
data for the post-cotunnite phase shown in Fig. [5] lends 
strong support to the viability of the P2i/m structure 
as a candidate for the post-cotunnite phase of Te02. We 
note that the impurity line indicated by a black arrow in 
Fig. [2] for the post-cotunnite phase is absent in the theo- 
retical data. Sato et al.^^ comment that other peaks in 
the experimental data may also be impurity lines, which 
could explain why some of the peaks are missing in the 
theoretical data. 

The CN of the Cmcm and Pl^jm structures are both 
found to be ten, in comparison with cotunnite which has 
CN=9. The average of the nine nearest-neighbour Te-0 
distances in the cotunnite structure at 130 GPa is 2.16 
A with a range of 2.02-2.28 A. For P2i/m the average 
is 2.23 A with a significantly larger distribution of bond 
lengths of 1.98-2.99 A. The Cmcm structure has Te-0 
bond lengths of 2.02-2.51 A with an average of 2.22 A, 
slightly less than for P2i/m due to the smaller distor- 
tions of the oxygen polyhedra in Cmcm. The oxygen 
coordination numbers with respect to neighbouring Te 
atoms are 4 and 5 in cotunnite and 4 and 6 in Cmcm 
and P2i/m. 

The P2i/m structure was studied in several other diox- 
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FIG. 5: (Color online) Band gaps of the cotunnite, P2i/m, 
Cmcm and PI structures as a function of pressure. 



ides to establish whether it might be a more general post- 
cotunnite phase. Ti02, P0O2, Th02, Se02, Si02, and 
Hf02 were tested, but no evidence was found to suggest 
that P2i/m is more stable than cotunnite in any of these 
materials. 



ELECTRONIC STRUCTURE OF THE PHASES 

The pressure dependence of the calculated band gaps 
of the structures are shown in Fig. [5l Above about 135 
GPa, the band gaps of the P2i/m, Cmcm, and PI struc- 
tures decrease with increasing pressure, however the co- 
tunnite band gap unexpectedly begins to increase sharply 
from a minimum of 0.49 eV before levelling off at higher 
pressures. This kink in the pressure dependence of the 
band gap of cotunnite approximately coincides with the 
change in compressibility seen in Fig. |31 The band gap 
of the Cmcm structure falls almost to zero by 250 GPa, 
although increasing the pressure further does not lead 
to overlapping valence and conduction bands. The insu- 
lating nature of the P2i/m phase is in agreement with 
the experimental observation of Sato et al.^ that the 
material was not opaque up to the highest experimental 
pressure of 150 GPa. 

The calculated band structures of the cotunnite and 
P2i/to phases at 130 GPa are shown in Fig.[6l The much 
larger band gap of the P2i/m phase is apparent. The 
eight lowest-energy bands for both phases arise from the 
02s states, and these lie at higher energies for the P2i/m 
structure, which has a somewhat smaller occupied band- 
width than cotunnite. Note that the band structure cal- 
culations were performed at the PBE-GGA level and are 
therefore expected to underestimate the true band gaps. 
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FIG. 6: (Color online) Band structures of the cotunnite (left) 
and P2\/m phases (right) at 130 GPa. The bands with ener- 
gies below zero are occupied and those at higher energies are 
unoccupied. 



CONCLUSIONS 

We have searched for the post-cotunnite phase of Te02 
using the AIRSS method. Our study supports the exper- 
imental observation of a post-cotunnite phase of Te02 
at pressures readily accessible within a diamond anvil 
cell. We predict a transition to a tenfold coordinated 
P2i/m phase at 130 GPa (at zero temperature), for 
which the calculated x-ray diffraction data are in good 
agreement with experiment. Although the P2i/m phase 
has a smaller volume than the cotunnite phase up to 
about 160 GPa, cotunnite has a slightly smaller volume 
at higher pressures, and we predict that the cotunnite 
phase re-enters at about 260 GPa. The P2i/m phase 
does not appear to be a general post-cotunnite phase for 
the dioxides. The P2i/m phase is found to be an in- 
sulator over the range of pressures studied, up to 400 
GPa, and hence should not appear opaque, in agreement 
with experiment, [sj Higher quality x-ray diffraction data 
are required to confirm whether our assignment of the 
P2i/m structure to the post-cotunnite phase of Te02 is 
correct. 
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